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ABSTRACT 

o' 

We present first resuhs from a narrowband imaging program for intermediate 
! redshift emission-hne galaxies using the newly commissioned FourStar infrared 

^ , camera at the 6.5m Magellan telescope. To enable prompt identification of Ha 

i-{ ', emitters, a pair of custom 1% filters, which sample low-airglow atmospheric win- 

^ ! dows at 1.19 /im and 2.10 /im, is used to detect both Ha and [OII]A3727 emission 

from the same redshift volume at 2=2.2. Initial observations are taken over a 
130 arcmin^ area in the CANDELS-COSMOS field. The exquisite image quality 
resulting from the combination of the instrument, telescope, and standard site 
conditions (~0'.'55 FWHM) allows the 1.19 fj,m and 2.10 /um data to probe 3a 
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emission-line depths down to 1.0x10"^^ erg s~^ cm~^ and 1.2x10"^^ erg s~^ cm~^ 
respectively, in less than 10 hours of integration time in each narrowband. For Hct 
at z = 0.8 and z = 2.2, these fluxes correspond to observed star formation rates of 
~0.3 and ~4 Mq yr~^, respectively. We flnd 122 sources with a 1.19 /xm excess, 
and 136 with a 2.10 fim excess, 41 of which show an excess in both bands. The 
dual narrowband technique, as implemented here, is estimated to identify >80% 
of 2; = 2.2 Ha emitters in the narrowband excess population. With the most 
secure such sample obtained to-date, we compute constraints on the faint-end 
slope of the z — 2.2 Ha luminosity function. Fitting of a pure power-law gives 
a — — 1.85±0.31, which is steeper than other recent estimates based on coarser 
selection techniques, but consistent within the typically large uncertainties that 
currently characterize such measurements. Combining our LF points with those 
at higher luminosities from other work, the slope decreases to a = —1.58 ±0.40. 
These "narrow-deep" FourStar observations have been obtained as part of the 
larger NewHa Survey, which will combine the data with "wide-shallow" imaging 
through a similar narrowband filter pair with NEWFIRM at the KPNO/CTIO 
4m telescopes, to enable study of both luminous (but rare) and faint emission-line 
galaxies in the intermediate redshift universe. 

Subject headings: galaxies: evolution — galaxies: photometry — infrared: galax- 
ies — surveys — techniques: photometric — galaxies: star formation 



1. Introduction 

The rest-frame optical emission lines are indispensable diagnostics of the activity in 
galaxies. In the local Universe, the measurement of lines from [Oil] A3727 through [SII]AA6716, 31 
has traditionally provided essential information on properties such as the star formation rate 
(SFR); the physical state of the ISM (e.g., electron temperatures and densities, gas-phase 
metallicities) ; the dust reddening; and may also indicate whether an AGN is present. Thus, 
it may seem like a conspiracy of nature that for 1 < 2; < 3, where star formation and AGN 
activity reach their global maxima (Shapley 2011 and references therein) and observing these 
lines is critical to the study of galaxy evolution, they are shifted into the near-infrared, where 
spectroscopic observations from the ground are difficult due to the bright sky background. 



^This paper is primarily based upon data gathered with the 6.5 meter Magellan Telescopes located at 
Las Campanas Observatory, Chile. 
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One way to partially circumvent the difficulties posed by the forest of near-infrared 
terrestrial lines is to use narrowband filters which are designed to transmit in dark regions of 
the sky spectrum. Such filters can be employed to blindly detect and measure emission lines 
at select redshifts. In the past, this technique has been successfully used with red-sensitive 
optical imagers to search for Lja emitters at the highest redshifts (e.g., Taniguchi et al. 2003 
and references therein). However, the near-IR instrumentation required to enable analogous 
surveys for galaxies via their rest-frame optical emission lines at "cosmic high noon" became 
readily available only around the mid-2000's. Initial attempts with the first generation of 
IR imagers with small fields-of-view were heroic, but had limited success. For example, 
Tephtz et al. (1998) collected a total of ~48 hours of Keck observations with NIRC (38" 
FoV), and found a sample of just 13 candidates. Even up to 2007, only a couple of hundred 
measurements of Ha in field galaxies past z~ 0.4 had been published (e.g., Tresse et al. 
2002; Erb ct al. 2006c). With the gradual maturation of wide- field near-infrared cameras 
however, efforts by a few independent groups (Villar et al. 2008; Geach ct al. 2008; Sobral 
et al. 2009; Hayes et al. 2010), including our own (Ly et al. 2011; J.C. Lcc ct al. 2012, in 
preparation), have now succeeded in increasing the total number of candidate Ha emitters 
at z > 0.8 by an order of magnitude. 

Nevertheless, while using strategically designed filters allows for deep samples of emission- 
hne galaxies at intermediate redshifts to be obtained, additional data are still needed to 
identify the line responsible for the detected narrowband photometric excess. Unambiguous 
determinations can be made with spectroscopy. However, since spectra are not usually im- 
mediately available for the selected sources, other strategies, typically based on broadband 
photometry and colors, are first used to deduce the nature of the excess (e.g., Fujita et al. 
2003; Kodama ct al. 2004; Ly et al. 2007; Shioya et al. 2008; Sobral et al. 2012). For 
surveys targeting intermediate redshift emission-line galaxies, such coarse filtering may be 
the only method available to discriminate between different line emitters, since a prohibitive 
investment of telescope time is still required to obtain near-lR spectroscopy for large and 
deep samples (e.g., Geach et al. 2008). The resulting contamination rate and incompleteness 
will vary depending on the depth of the photometry and number of colors available (e.g., 
15-30% contamination if Riz selection is used to select Ha emitters a,t z — 0.8; Ly et al. 
2011), but will be most severe for young, blue galaxies, since the primary spectral feature 
providing information on the approximate redshift, the Balmer/4000 A break, is weak or 
non-existent in such populations. 

To address such issues, we have attempted to improve upon the basic narrowband 
technique. Specifically, we wish to enable the prompt identification of Ha emitting galaxies 
at z ~ 2, and also provide an initial glimpse into the state of the ISM which is encoded into 
the rest-frame optical emission-lines. Rather than use a single narrowband filter, we have 
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designed a pair of filters, both of which sample dark regions of the near-IR sky spectrum, 
and are coupled to capture Ha and [OII]A3727 over the same redshift volume. Emission- line 
galaxy candidates are selected through detection of a photometric excess in a filter at 2.10 
/im, and Ha emitters are then isolated through the presence of a second narrowband excess 
at 1.19 /im|^ The Ha line is of particular interest because its fiux is a primary indicator of 
the instantaneous star formation rate (SFR) - the dust-corrected Ha luminosity is directly 
proportional to the ionizing fiux. The [OH] fiux has also been used as an SFR indicator, 
but its strength is more strongly dependent on the metallicity and the ionization state of 
the gasjfl The observed [OH] /Ha ratio thus can provide a coarse combined measure of the 
metallicity and reddening (e.g., Nakajima et al. 2012). Furthermore, this technique yields an 
excellent sample of emission-line galaxies at z=2.2 for follow-up NIR spectroscopy, as three 
key lines, [OH], and Ha along with the fianking [Nil] lines, are already pre-selected to fall 
in clean windows of the OH sky spectrum. 

Here, we describe the narrowband imaging in this filter pair obtained by the NewHa Sur- 
vey (J.C Lee et al. 2012, in preparation) using the newly commissioned FourStar infrared 
camera (S.E. Persson et al. 2012, in preparation) at the 6.5 Magellan-Baade telescope. 
Initial observations are carried out over an ll'xll' area in the Cosmic Evolution Survey 
field (COSMOS; Scoville et al. 2007) which overlaps HST WFC3/IR and ACS coverage by 
the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS; Grogin et 
al. 2011; Koekemoer et al. 2011), and WFC3/IR grism observations by the 3D-HST Survey 
(van Dokkum et al. 2011). We take advantage of the wealth of deep, multi- wavelength imag- 
ing observations available in COSMOS to construct spectral energy distributions (SEDs), 
and compute photometric redshifts for the narrowband excess sources. A unique aspect 
of the ensemble of data used in the photometric redshift analysis is that it includes near- 
IR medium-band imaging through five bands from the FourStar Galaxy Evolution Survey 
(Z-FOURGE; I. Labbe et al. 2012, in preparation). The Z-FOURGE filter set provides a 
powerful complement to the NewHa filter pair, as it was designed to locate the Balmer- or 
4000 A- breaks in galaxies at 1.5 < 2 < 3.5, a range in which a significant fraction of 2.10 /im 
narrowband selected emission-line galaxies are expected to lie. Together with other public 
datasets, photometry in 28 filters is used, spanning from the U to Ks bands (^3800 A to 
~2.2 /im). With the SEDs and photometric redshifts, we examine the composition of the 
narrowband excess samples, and test the effectiveness of the dual narrowband technique. 



^This [Oil] /Ha dual narrowband selection technique has just recently been also used by Sobral et al. 
(2012) to obtain a galaxy sample at z = 1.47. There, [Oil] and Ha are detected using Subaru/Suprime-Cam 
imaging in NB921, and UKIRT/WFCAM imaging in "iVB//" (A=1.617 /xm), respectively. 

•^See Kennicutt (1998) and Kennicutt & Evans (2012) for a review of star formation rate indicators. 
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The remainder of the paper is organized as follows. In Section 2, an overview of the 
NewHct and Z-FOURGE parent surveys is given, and the FourStar observations and data 
reduction is summarized. The method for selecting narrowband excess samples and identify- 
ing dual narrowband emitters is described in Section 3. Section 4 discusses the construction 
of SEDs, and the computation of photometric redshifts. An inventory is taken of the various 
detected line-emitters using the photometric redshifts and BzK colors in Section 5, and the 
effectiveness of the dual narrowband excess method for identifying a clean sample of 2; = 2.2 
Ho; emitters is evaluated. In Section 6, this sample is compared with the rest-frame UV- 
selected z ~2 galaxies of Erb et al. (2006a,b,c) to (i) illustrate the parameter space probed 
by our observations relative to a well-studied sample of intermediate redshift star-forming 
galaxies; and (m) test the applicability of local scaling relations used to remove the contri- 
bution of [NII]AA6548,83 to the narrowband excess flux, and correct for dust attenuation. 
Initial constraints on the z = 2.2 Ha luminosity function based on our new observations 
are also presented. Finally, we summarize our findings, and describe future work with the 
dataset in Section 7. 

Throughout the paper, the standard ACDM cosmology with [Q/^, Qm, ^70] = [0-7, 0.3, 
1.0] is assumed. Magnitudes are reported on the AB system. Star formation rates are 
computed using the Kennicutt (1998) prescription and hence assume a Salpeter IMF with 
mass hmits of 0.1 and 100 Mq. 

2. Observations and Data Reduction 

2.1. The New Ha Survey 

The FourStar narrowband imaging used in this analysis has been obtained as part of the 
overall NewHo; Survey, a program which extends searches for emission-line galaxies into the 
intermediate redshift universe (J.C. Lee et al. 2012, in preparation). The NewHo; Survey has 
been designed to capture statistical samples of both luminous (but rare), and faint emission- 
line galaxies, by combining the near-infrared imaging capabilities of NEWFIRM (Probst et 
al. 2008) at the KPNO/CTIO 4-m telescopes (FoV 27'.6 x 27f6) to cover large areas, and the 
exquisite image quality (~0'.'5-0'.'6 FWHM) of FourStar at the Las Campanas Magellan 6.5-m 
(FoV 10^8 X lO'.S) to probe luminosities that are a factor of three deeper over smaller fields. 
Follow-up programs with GALEX, Spitzer and the Magellan-IMACS spectrograph are in 
various stages of progress as described in (J.C. Lee et al. 2012, in preparation) Narrowband 
observations at 1.19 //m with NEWFIRM are complete, and results on (i) the Ha luminosity 
function and SFR volume density at ^ = 0.8, and (ii) the dust attenuation oi z — 0.8 
Ha emitters based on Balmer decrements measured from IMACS optical spectroscopy are 
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reported in Ly et al. (2011) and Momcheva et al. (2012), respectively. 

For both NEWFIRM and FourStar, a pair of custom narrowband filters, which sample 
low OH-airglow windows in the near-IR sky spectrum, has been manufactured. The band- 
passes are made to fit within high atmospheric transmission regions at 1.19 and 2.10 fim 
which are relatively free of terrestrial sky lines over 1% in A (Figure [1], top). The filters 
are sensitive to Ha at z=0.8 (near the beginning of the ten-fold decline in the cosmic SFR 
density) and at z=2.2 (within the peak of the cosmic star formation history). Fortuitously, 
filters centered around these two wavelengths can also be designed so that both [OII]A3727 
and Ha emission can be detected from the same galaxies at z = 2.2 (Figure [2]). Again, such 
dual narrowband excess detections help to differentiate Ha emitters from other emission-line 
sources prior to confirmation with near-IR spectroscopy. Another useful feature of the par- 
ticular redshift singled out by this filter pair is that Lya appears at ~3900 A, which is just 
above the blue cut-off of most optical cameras. Thus, a third narrowband can be designed 
to also detect Lya in the same volumes, as has been done by Nakajima et al. (2012) (see 
also Hayes et al. 2010). 

The central wavelength and FWHM of the two filters built for FourStar are 11908 A 
and 137 A and 20987 A and 208 A. Hereafter, we refer to these filters as NB119 and NB210, 
respectively. 

The observational goal of the NewHa survey's FourStar program is to probe line lumi- 
nosities that are roughly an order of magnitude fainter than the "knee" of the Ha luminosity 
function at z = 2.2, to allow for example, robust constraints on the faint end slope to be 
derived. Based on measurements made by Geach et al. (2008) and Hayes et al. (2010), this 
depth is estimated to be around Ha luminosities which translate into to observed SFRs of 
5 M0 yr~^ at z ~ 2.2, (fiuxes of ~ 1.7 x 10"^''^ ergs s~^ cm~^). This limit can be reached 
in NB210 imaging with a limiting AB magnitude of ~23.7. The corresponding depth re- 
quired in the NB119 imaging to detect [Oil] for the faintest Ha z = 2.2 emitters is roughly 
8.5 X 10"^^ ergs s~^ cm~^ (25.2 AB), if it is assumed that the observed Ha/[OII] ratio is 2. 
The assumption is based on spectroscopy of local star-forming galaxies. Such studies show 
an average ratio of ~2, with values ranging between 1 and 4 for massive L* systems, and an 
average of ~1.5 with values between 1 and 2 for lower luminosity dwarf galaxies (Jansen et 
al. 2001; Moustakas et al. 2006; see also FigureE]). The effectiveness of these chosen relative 
depths for identifying the 2; = 2.2 Ha emitters in the NB210 excess sample is examined in 
Section [51 
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2.2. Z-FOURGE 

Z-FOURGE, the FourStar Galaxy Evolution Survey (I. Labbe et al. 2012, in prepa- 
ration), is a program which aims to characterize the properties of galaxies at ^ > 1 by 
performing the first crucial (and often limiting) step of deriving reasonably accurate red- 
shifts (i.e., -2% in ^). Ultimately, Z-FOURGE seeks to measure such redshifts for a 
sample of ~30,000 i^s-selected galaxies at 1-5 < 2 < 3.5. The principal focus of Z-FOURGE 
is to investigate the red, early-type population using mass-selected samples. Since quiescent, 
non-active L* galaxies at 2; > 1 are too faint for spectroscopy, Z-FOURGE's strategy is to 
use a custom set of five medium-band, near-infrared filters to constrain the location of the 
Balmer/4000 A break. The filters sample wavelengths between 1 to 1.8 /im, and splits the 
standard J and H bands into three (Jl, J2, J3) and two {HI, Hs) bands, respectively. 
A similar filter set has also been used with NEWFIRM by the NEWFIRM Medium-Band 
Survey (van Dokkum et al 2009; Whitaker et al. 2011). 

Bandpasses for the Z-FOURGE medium-band filter set are shown together with the 
NewHa filter pair in Figure [1] 

2.3. FourStar Observations 

Pilot study observations for the NewHa and Z-FOURGE surveys were carried out during 
"science demonstration time" in FourStar's first semester (2011 A) on the Magellan-Baade 
6.5m telescope. Both programs targeted the same region centered on a = 10^00™'32'^.4 
5 = +02°16'58", within the CANDELS-COSMOS field. Standard near-infrared observing 
procedures were followed. The narrowband observations are summarized here. A description 
of the medium-band observations can be found in the Z-FOURGE survey paper (I. Labbe 
et al. 2012, in preparation). 

Narrowband observations were performed over a span of eight nights (24 January; 19, 
21-22 February; 13-14, 19-20 March 2011), which were shared with the Z-FOURGE survey. 
Cumulative integration times of 9.65 and 8.39 hrs in the NB119 and NB210 filters were 
obtained on a single FourStar pointing. The exposure times for individual NB119 and NB210 
frames were 4 and 2 minutes, respectively. The data were taken with the instrumental gain in 
its low noise mode (~1.3 e-/ADU). Semi-random dither patterns, where the spacing between 
adjacent frames is maximized and the pointing centers are distributed over a 27" box, are 
used. The dithers serve to bridge the 18" gaps between FourStar's four 2048 x 2048 HAWAII- 
2RG detectors; enable filtering and removal of bad pixels; and perhaps most importantly, 
provide a measurement of the sky in each pixel for accurate sky-subtraction. The median 
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seeing during the observations was 0'.'6, and varied between 0'.'45 and O'.'QO. Hence, the PSF 
was adequately sampled by FourStar's 0'.'159 pixels. Sky conditions were mostly photometric, 
although some data were taken through thin cirrus. Continuum measurements are provided 
by imaging in the Kg band for the NB210 data, and by the J2 and J3 bands for NB119. 

A summary of the observations in given in Table [H 

2.4. FourStar Data Reduction 

The FourStar data were reduced using a custom pipehne, adapted from the IDL code 
used for the NEWFIRM Medium-Band Survey (see Whitaker et al. 2011). The reduction 
follows standard, iterative (two-pass), near-infrared reduction techniques (e.g., Labbe et al. 
2003) to subtract the sky background; reject artifacts; optimally weight and stack the indi- 
vidual frames; and is consistent with the procedures used to process the NewHa narrowband 
data taken with NEWFIRM (Ly et al. 2011; J.C. Lee et al. 2012, in preparation). Cor- 
rections are applied for non-linearity and geometric distortion, and flat-fielding is currently 
performed using twilight flats taken in each band. Comparison with independent, publicly 
available, broadband near-IR imaging in COSMOS shows that large-scale systematic varia- 
tions in the sensitivity over the array have been reduced to <5%. Photometric calibration 
is performed using observations of the spectrophotometric standards GD71 and GD153, and 
the resultant zeropoints are accurate to within ~2%. Further optimization of the reduction 
and calibration is in progress. 

The 3-0" limiting depths of the stacked imaging in a 1'.'2 diameter aperture are given 
in Table [T] A 1'.'2 diameter aperture is roughly twice the size of the PSF, and will contain 
99% of light from a point source, assuming a Gaussian profile. It is the size of the aperture 
used to select narrowband excess sources, as described in the next section. The narrowband 
observations successfully achieve the NB210 target depth described in Section [2?T] (1.7x 10~^^ 
ergs s~^ cm~^) at the 5-cr level, and surpass it at the 3-0" level. The NB119 observations 
reach the target depth at about the 3-0" level. Observations in the continuum bands reach 
flux limits which are approximately twice as deep. The limits are computed by measuring 
the background in a large number of random apertures. 
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3. Sample Selection 

3.1. Selection of NB210 Excess Objects 

Sources that show a significant i^^— NB210 color excess are selected as emission-line 
galaxy candidates. The selection procedure follows general techniques commonly used in 
narrowband surveys (e.g., Villar et al. 2008, Sobral. et al. 2009, Ly et al. 2011, Nakajima 
et al. 2011). 

First, we detect sources in the NB210 image with SExtractor, using the default param- 
eter values which require that detections have 5 adjacent pixels which are each individually 
above 1.2cr of the background. Exposure maps are used as weight images to avoid spurious 
detections in low signal-to-noise regions. A Gaussian convolution filter with a FWHM of &!6 
is applied to improve identification of faint, extended, objects. A total of 4865 sources are 
found. Ks photometry is measured for each of these NB210 detections in matched apertures 
using SExtractor in "dual-image" mode. As discussed in the previous section, an aperture 
of 1'.'2 diameter is chosen to compute the color. 

Using this photometry, NB210 excess sources are then selected by requiring that: 

(i) A(/Cs-NB210), the Ks-NB210 color corrected for the median color of bright continuum 
objects with -0.4 < irs-NB210< 0.4 and 19.0 <NB210< 21.5, has a positive value which is 
significant at the 3a level or higher (curve in Figure E]); and 

(ii) A(i^'^-NB210) > 0.275, which is five times the scatter in A(ir^-NB210) of the sample 
of continuum objects defined in (i) (upper horizontal line in Figure [3]). The minimum color 
represents a threshold in the rest-frame equivalent width of Ha at ^ = 2.2 of 20 A (cf. Ly 
et al. 2011, eq. 11). 

The selection results in a sample of 139 objects, 13 of which are undetected in the Kg- 
band. The objects are individually inspected on both the narrow and broad-band FourStar 
images, and none appear to be spurious detections. SEDs constructed from broadband 
photometry (Section H]) are examined to identify stars. Three stars are found and removed, 
reducing the NB210 excess sample to 136 objects. 

3.2. Selection of Dual Narrowband Excess Emitters 

To determine which NB210 excess sources are "dual emitters" (i.e., they also show a 
photometric excess in the NB119 filter, and are likely to be Ha emitters at z = 2.2), we 
perform photometry on the NB119 image, and on a combined J2 and J3 (hereafter "J2J3") 
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image. The measurements are again made with SExtractor in dual-image mode, with the 
detection performed on the NB210 image. This produces catalogs of NB119 and J2J3 
photometry at the locations of each of the 4865 NB210 detections. 

Criteria analogous to those described in the previous section are used to select the NB119 
excess sources from these catalogs. The corresponding color-magnitude diagram is given in 
the right panel of Figure [31 The fraction of the NB210 excess sample which also show an 
excess in NB119 is 30% (N=41). In both panels of Figure [3], the NB119 excess sources are 
enclosed with blue circles, and the NB210 excess sources with red circles; hence the points 
marked in both red and blue are dual emitters. All of the dual emitters are detected in the 

3.3. Selection of NB119 Excess Objects 

Following the above procedures source detection on the NB119 image and matched 
aperture photometry on the J2J3 image can also be performed to obtain an NB119 excess 
selected sample without reference to the NB210 data. A total of 6894 sources are detected, 
with 122 showing A( J2 J3-NB119)>0.21 at at least 3a significance (Figure S]). 

4. SEDs and Photometric Redshifts 

SEDs for the narrowband detected sources are constructed from the ensemble of publicly 
available optical broad-band and medium-band data in the COSMOS field, in combination 
with the Magellan-Fourstar NewHa and Z-FOURGE imaging described above. The SEDs 
are fit using the EAZY code (vl.1.10; Brammer et al. 2008) to derive photometric redshifts. 
We use the same optical dataset and a procedure similar to those in Whitaker et al. (2011), 
who performed a near-IR medium-band study of the COSMOS field. With EAZY, linear 
combination of seven template spectra, which span a broad range of galaxy ages, along with 
a Bayesian prior on the observed total Kg magnitude, are used to to fit the observed SEDs. 
Based on these fits, a variety of different characteristic photometric redshifts are calculated 
from the probability distribution of redshifts for a given galaxy. In the analysis that fol- 
lows, ^^Zpeak" the redshift peak in the distribution with the highest integrated probability, is 
adopted as the photometric redshift. 

The public datasets in the COSMOS field that are included in this analysis consist of 
imagmg m 22 filters from the U to Z bands (~3800 A to 9000 A). We use data from: (1) the 
Deep Canada-France-Hawaii Telescope Legacy Survey (CFHTLS ugriz, Erben et al. 2009; 
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Hildebrandt et al. 2009), and (2) Subaru Suprime-Cam {BVr'i'z', Taniguchi et al. 2007; and 
12 medium-band filters from 4200 A to 8200 A, Y. Taniguclii et al. 2012, in preparation). 

Combining tfie COSMOS optical imaging in 22 filters with FourStar near-IR imaging 
in 5 Z-FOURGE medium-bands, and a standard Kg filter, provides up to 28 points for each 
individual SED. The imaging in the various bands were all convolved to the same PSF (1'.'24 
FWHM), and photometry was measured in I'.'S diameter apertures. Details on the PSF 
matching of the imaging, and the generation of a matched multi-wavelength flux catalog are 
given in Spitler et al. (2012) and I. Labbe et al. (2012, in preparation). 

Our photometric redshifts are in reasonable agreement with available spectroscopic red- 
shifts for the NB210 excess sources. We searched the zCOSMOS catalog (Lilly et al. 2007) 
for spectroscopic redshifts, and found 4 NB210 excess emitters with matches within 1". All 
had Zspec < 1; and none were dual excess emitters. Each of the 4 galaxies have individual 
photometric and spectroscopic redshifts which are consistent to within 1% in -r-^ — . The 
NB210 excess for two galaxies can be identified as originating from Pact aX z = 0.12, and 
from HeIA10830 at z = 0.94 for the other two. 

An additional check of the accuracy is also performed by comparing the photo-z's of 
the larger parent sample of NB210 detections (as opposed to the NB210 excess detections) 
with available spectroscopic redshifts. There are 458 galaxies with NB210 detections that 
have reliable spectroscopic redshifts (zCOSMOS "confidence class" > 90%) within 1". Com- 
parison of the photometric and spectroscopic redshifts for Zgpec < 1.4 show a dispersion of 
(^ij^ — )=0.034, with no significant mean offset. This accuracy may appear to be low rela- 
tive to the ~1% reported by Ilbert et al. (2009) for their own set of photo-z's based on 30 
bands of photometric data. However, the NB210 detections which have spec-z's available in 
zCOSMOS are much fainter on average (the median in a 3" aperture is 22 mag) than the 
17.5 < < 22.5 sample for which the 1% photo-z's have been derived. Photo-z's computed 
by Ilbert et al. (2009) are available for 90% of the NB210 detections with spec-z's, and when 
these are compared, cr(—^ — )=0.041, consistent with the accuracy of the photo-z's derived 
here. 



5. Composition of the Narrowband Excess Samples 

5.1. NB119 Excess Emitters 

The photometric redshifts provide a first look at the composition of the narrowband 
excess samples. As expected, the photo-z distribution shows a series of peaks at redshifts 
where strong emission-lines appear in the filter. Figure [5] shows that the main emission-lines 
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detected in the FourStar NB119 excess sample (N=122) are: 

• [SIII]AA9069,9532 at z = 0.31 and 0.24, 

• [SII]AA6717,6731 at z = 0.76, 

• Ha at z = 0.80, 

• [OIII]AA4959,5007 at ^ = 1.39 and 1.36, 

• at 2; = 1.44, and 

• [OII]A3727 at 2 = 2.19. 

The majority of the sample is at 2; > 1. Roughly a quarter are [OIII] or emitters, 
while a third are [Oil] . The distribution is considerably different from, but highly comple- 
mentary to those of previous narrowband surveys at 1.19 /xm. Past studies reached half of 
the depth of the current observations or less, but covered larger fields, and hence were dom- 
inated by foregroimd Ha emitters at z ~ 0.8 (e.g., Villar et al. 2008; Sobral et al. 2009; Ly 
et al. 2011). Here, Ha or [Sll] candidates comprise only about 10% of the sample (N ~ 15), 
but nearly all of these objects have luminosities lower than 10^"^ ergs s^'^, consistent within 
the expectations based on extrapolations of published Ha luminosity functions a.t z = 0.8, 
and the volume covered in this one FourStar pointing. An Ha luminosity of 10^^ ergs s""^ 
corresponds to 0.8 Mq yr~^ , and the candidates span SFRs down to ~0.3 Mq yr~^ , similar 
to the present activity of the Large Magellanic Cloud (Whitney et al. 2008; Harris & Zarit- 
sky 2009). Ultimately, the NewHa FourStar NB119 data will allow us to derive the most 
stringent constraints on the faint-end slope of the Ha luminosity function, and probe the 
star formation properties of true dwarf galaxies at z — 0.8. 

In principle, detection of Lya at z — 8.8 with the NB119 filter is also possible, but 
unlikely given the expected number densities of Lya emitting galaxies at early times, and 
the volume and depth probed by our data (cf., Nilsson et al. 2007; Sobral et al. 2009). If 
the Lya luminosity function at z = 6.6 (Ouchi et al. 2010) is assumed to hold at z = 8.8, 0.4 
Lya emitters would be expected in one FourStar NB119 pointing with depth of 10^^^ ergs 
s^^ cm^^. However, this assumption is likely to lead to an overly optimistic estimate as Ono 
et al. (2012) has reported that the Lya luminosity function appears to decrease by a factor 
of two to three from 2; ~ 6 to 7. We are currently investigating whether any of the NB119 
sources are Lya candidates, and results will be reported in a forthcoming paper. Further 
discussion of the NB119 excess selected sample is deferred to future work. The remainder of 
this paper will focus on the NB210 excess sample, and dual NB210/NB119 excess emitters. 
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5.2. NB210 and Dual Narrowband Excess Emitters 

In comparison, ttie NB210 filter is sensitive to a broader array of emission-lines. The 
filter will not only detect lines redward of 1.19 /um (such as Paa, Pa/3 and all the high order 
Brackett lines in between), but also those in the near-IR below 1.19 fim (such as HeIA10830 
and Pa7) since the NB210 filter will capture these lines at higher redshift and thus probe 
significantly larger volumes relative to the NB119 filter. As shown in Figure |6l the main 
emission-lines detected in the NB210 excess sample (N=139) are: 

• Brackett lines from Bid to the series limit from z = 0.08 to 0.43, 

• Paa to Pae at ^ = 0.12, 0.64, 0.92, 1.09, 1.20 (series limit at z = 1.56), 

• [FeII]A16435 at ;z = 0.28 

• Hel A10830 at ^ = 0.94, 

• [S111]AA9069,9532 at z = 1.31 and 1.20, 

• [SII]AA6717,6731 at z = 2.12, 

• Ha at 2; = 2.19, 

• [Olll]AA4959,5007 at ^ = 3.23 and 3.19, 

• H/3 at z = 3.31, and 

• [OII]A3727 at z = 4.62, 

The majority of NB210 excess emitters are also at 2; > 1. Over a third are Ha or [SII], 
about a quarter are [OIII] or H/3, and there may possibly be a few [OH] emitters in the 
sample. Approximately 40% are due to H or He recombination emission from galaxies in the 
foreground. There is also a spike at z ~0.3 containing about 15 objects, where the shock 
excited [Fe II] line is expected to appear in the NB210 filter. As mentioned above, only four 
objects have spectroscopic redshifts from the zCOSMOS survey. These show that Paa at 
z = 0.12 is responsible for the NB210 excess in two galaxies, and that the excess is from 
HeIA10830 at ^ = 0.94 in the other two. 

The photo-z distribution of NB210 excess emitters that also have an excess in the 
NB119 band is shown with the shaded histogram in Figure |6l About 90% (36/41) of the 
dual narrowband excess emitters have photo-z's between 1.8 and 2.6, demonstrating that the 
technique succeeds in identifying the z = 2.2 Ha emitters in the sample. Five dual excess 
emitters have photo-z's less than 0.4, but examination of their SEDs and photo-z probability 
distribution functions, shows that they can also be reasonably fit with models at z = 2.2. 
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The effectiveness of tlie NB210/NB119 dual narrowband excess teclinique for selecting 
a clean sample of Ha at z = 2.2 can also be shown by plotting the objects on a BzK color- 
color diagram. The BzK diagram has frequently been used to obtain samples of galaxies with 
redshifts between about 1.4 and 2.5 with relatively high completeness and low contamination, 
and also to distinguish between star-forming and passively evolving galaxies, as well as 
stars (e.g., Daddi et al. 2004; Kong et al. 2006; McCracken et al. 2010). The method 
was introduced, calibrated and tested by Daddi et al. (2004) using both spectroscopic 
samples and stellar population synthesis models. In Figure [71 the original criteria given 
by Daddi et al. (2004) are marked with solid lines0 All but three of the dual excess 
emitters (blue solid symbols) satisfy the criteria for classification as intermediate redshift 
star-forming galaxies. Even the three objects that fail to meet the criteria fall in a region 
slightly below the star-forming boundary, {z — K) — {B — z) > —0.2, in which populations 
dominated by very young populations (< 10^ yrs) can exist. Daddi et al. (2004) notes that 
to capture such galaxies, the boundary must be extended to {z — K) — {B — z) > —0.8 
(dashed line), but that this comes at the cost of higher contamination by populations at 
z < 1.4. The five dual emitters that have photo-z's less than 0.4 (triangles) straddle the 
formal {z — K) ~ [B — z) > —0.2 star-forming BzK boundary, with two falling outside the 
boundary. Optical spectroscopy can be performed to check whether these are in fact low 
redshift galaxies. In any case, the contamination of the dual excess sample with non-Ha 
emitters will be very low, < 5%. 

FigureElalso shows that there are NB210 excess objects that have photo-z's which would 
make them Ha candidates, but do not have NB119 excess detections. If the photo-z range 
for Ha is taken to be between 1.8 and 2.6 (as defined by the range of the majority of dual 
excess emitters), there are nine such objects. The two most likely explanations for these 
objects is that the NB210 excess is due to [SII]AA6717, 31, or that the depth of the NB119 
image was not sufficient to detect [OH]. 

To examine the latter possibility, the Ha+[NII] and [OH] fluxes for the objects in the 
NB210 excess sample are plotted in Figure El and compared with those observed locally from 
the Sloan Digital Sky Survey (SDSS) Data Release 7 sample (DR7; Abazajian et al. 2009). 
Emission-line fluxes are computed from the FourStar narrowband and broadband photometry 
in the usual way (e.g., see equations 9 and 10 in Ly et al. 2011). Spectroscopic emission-line 
fluxes for the SDSS sample are made publicly available by the MPA/ JHU group0. The fluxes 



''The empirically derived transformation given by McCracken et al. (2010) to account for differences 
between the Subaru Bj filter used to image the COSMOS field, and the B-VLT filter used by Daddi et al. 
(2004) is applied to the data shown. 

http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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have been extracted from Gaussian fits to the fine profiles, where the continuum has been 
subtracted using a best-fit model generated by a combination of simple stellar population 
at various ages from the Bruzual & Chariot spectral synthesis code (Bruzual & Chariot 
2003). Corrections for Galactic foreground reddening are applied. SDSS spectra for a total 
of 818,333 galaxies at z < 0.7 are available, but here we only consider the sub-sample of 
~60,000 star-forming galaxies with 0.005 < z < 0.22 (following Brinchmann et al. 2004) 
and an empirically estimated fiber aperture correction less than 2. The 3a sensitivities of 
the narrowband imaging to these lines is also illustrated by indicating which of the SDSS 
galaxies would be detectable if they were instead ai z = 2.2. Clearly, if the full range of local 
galaxies exists at redshift of two, current near-IR narrowband surveys, of which NewHa is 
one of the most sensitive to faint line emission, would still only be able to probe the most 
luminous Ha and [Oil] emitters. 

In Figure El all of the dual narrowband excess emitters lie in the region above both the 
NB119 and NB210 3a line fiux sensitivity limits as expected (yellow region). The region 
where objects would be only detected with a NB210 excess is indicated in cyan. Here, more 
uncertain NB119 excess fiuxes (which would correspond to [OH] if the objects were indeed at 
z = 2.2) down to the la limit (marked by the row of upper-limit symbols) are also plotted. 
It is clear that even without additional photo-z or SED information, some objects can be 
eliminated as Ha candidates because they would have [OH] /Ha ratios that are unphysical. 
However, the nine galaxies that do not have 3a excess detections in NB119, but yet have 
photo-z's which would make them Ha candidates, all would have [OH] /Ha ratios that fall 
within the range of physically plausible values. Further, all except one have nominal NB119 
excess detections above la, which would not be expected if the NB210 excesses were mostly 
due to [SII]. This suggests that the majority of these objects are Ha emitters at z = 2.2. 
Therefore, the dual narrowband technique, as implemented here, is estimated to identify 
>80% of z = 2.2 Ha emitters detected in the NB210 excess sample. Near-IR spectroscopy 
will be required to confirm this success rate. 



6. Constraints on the Ha Luminosity Function at z = 2.2 

With a clean, complete sample of Ha emitters in hand, constraints on the z = 2.2 Ha 
luminosity function can be computed. We follow the procedures discussed in detail in Ly et 
al. (2011), which presents our results on the Ha luminosity function at z = 0.8, using NB119 
imaging taken with the much wider FoV NEWFIRM camera on the KPNO 4m telescope. 
Here, we outline the procedure, pointing out differences between the analyses when they 
arise. 
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The key steps include (1) transforming the narrowband photometric excesses into dust- 
corrected Ha luminosities; (2) computing the effective volume of the observations; (3) es- 
timating and applying completeness corrections; and (4) fitting a function, usually that of 
Schechter (1976) (but see Salim et al. 2012, in preparation), to model the distribution. 
The step where methods which have been previously applied for analogous studies at lower 
redshifts cannot be directly transplanted for the current analysis is the first one. 

To compute the intrinsic Ha luminosity from the narrowband excess, two main cor- 
rections are required, one for emission from [NH]AA6548,83 (since the flux from these lines 
will contribute to the narrowband excess as R~100 filters are generally used), and another 
for dust attenuation within the galaxy. Because measurements of [NH]/Ha and A(Ha)|^ are 
not immediately available for the individual members of a sample, coarse empirical scaling 
relationships calibrated on local samples such as the SDSS have been employed to estimate 
their values (e.g., Villar et al. 2008; Sobral et al. 2009; Lee et al. 2009; Ly et al. 2011). 
For example, in previous Ha narrowband studies, relations providing [NII]A6584/Ha as a 
function of the rest-frame Ha equivalent width (Villar et al. 2008), and the intrinsic Ha 
luminosity as a function of the observed Ha luminosity (Hopkins et al. 2001) have been 
used. However, it is well known that such relationships have scatters of a factor of two or 
more (e.g., Villar et al. 2008; Kennicutt et al. 2008; Lee et al. 2009; Momcheva et al. 2012). 
Moreover, it is unclear whether these scaling relations will be valid at z = 2, as the gas 
metallicities and dust reddenings of galaxies of a given stellar mass are known to decrease 
with increasing redshift (Erb et al. 2006a; Reddy et al. 2006; Maiolino et al. 2008; Reddy 
2010). 

Though we do not yet have independent measurements of [Nil] /Ha and A(Ha) for 
our z = 2.2 Ha sample, we can take a first look at where these values are likely to lie via 
comparison to the well-studied rest-frame UV-selected z ~2 sample of Erb et al. (2006a,b,c), 
for which Ha spectroscopy has been obtained, and dust reddenings have been computed both 
from SED fitting and the UV continuum slope (Erb et al. 2006c). Such an exercise also serves 
to place the current dataset in the more familiar context of a broadly referenced sample of 
star-forming galaxies at z ~2. 

In FigureO the rest-frame equivalent width, and the line luminosity, scaled by 7.9x10"^^ 
to give approximate star formation rates (Kennicutt 1998), are plotted against the Ks mag- 
nitude for both samples. For the UV-selected sample, the measurements are made from slit 
spectroscopy, and so are for only the Ha line. These measurements include a factor of two 
correction to account for slit losses, as determined by Erb et al. (2006c). For the FourStar 



^attenuation of Ha emission in magnitude units 
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Ha emitters, the plotted values include contributions from both [Nil] and Ha. No dust 
corrections have been applied to either sample. Kg magnitudes on the AB scale are shown 
for both, where Kg [KS\=Ks [Vega]+1.83 is applied to the Erb data. 

Although the EW range (and hence range of specific SFR) spanned by the two samples 
are similar, it is clear that the Ha emitters probe SFRs and continuum luminosities (and 
hence stellar masses) that are significantly lower than the UV selected sample. The medians 
of the distributions for both quantities are at least a factor of two smaller for the Ha emitters. 
Hence, if the [Nil] /Ha values and reddenings derived for the UV-selected sample were to be 
used to formulate corrections for the Ha emitters, it seems that some extrapolation of the 
corrections would be required. 

We first examine the reddenings for the UV-selected sample found by Erb et al. (2006b) 
in the left panel of Figure [TDi Here, A (Ha) is plotted as a function of the observed SFR, 
to enable comparison with the locally calibrated Hopkins et al. (2001) relation. Erb et al. 
(2006b) provide E(B-V)'s derived by SED fitting, and these are scaled using the Calzetti 
et al. (2000) extinction law to give the attenuation in Ha. A(Ha) is computed both by 
assuming that the reddening in the gas is about twice that of the stars, which is found locally 
(Calzetti 2000); and, that it the same as that of the stars, which seems to be supported by 
intermediate redshift observational studies (Erb et al. 2006b, but see Forster-Schreiber et 
al. 2009). There is no correlation with the observed SFR, and in both cases, the Hopkins 
prescription over-predicts the attenuation. The median Ha attenuation is 0.5 mag when 
E(B-V)stars=E(B-V)gas is assumed, and so for the current analysis, we simply adopt this 
fixed value to perform the dust correction. 

Next, in the right panel of Figure [TOl the equivalent widths from Erb et al. (2006c) 
are used used to predict [NII]A6584/Ha using the local relation of Villar et al. (2008). 
The estimates are plotted as a function of Kg magnitude. This enables comparison with 
the Erb et al. (2006a) measurements, which were performed on long-slit spectra, stacked 
in bins of stellar mass. Kg magnitudes are provided for each mass bin. Figure [10] shows 
that the local relation actually yields reasonable predictions, albeit with large scatter. The 
consistency between the locally-based estimates and the z ^ 2 measurements may initially 
seem surprising given the known evolution of the mass-metallicity relation, but is likely a 
consequence of the "Fundamental Metallicity Relation," or FMR, proposed by Mannucci et 
al (2010). The FMR is a more general relation between mass, SFR and metallicity, with 
a small residual dispersion in the metallicity of only ~12% based upon the SDSS, and is 
shown to hold for galaxies to z ~2.5. Collapse in one dimension yields a correlation between 
metallicity and specific SFR, which is traced by the Ha EW. Here, although the Villar et 
al. (2008) scahng relation appears to be valid, we choose to use a linear fit to the Erb et 
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al. (2006c) measurements to estimate the [NII]A6584/Ha values for the Ha emitters from 
their K, magnitudes. The fit is given by [NII]A6584/Ha = -0.15K, + 3.67. For > 23.3, 
[NII]A6584/Ha = 0.04 is assumed. The vahdity of both this relation and the 0.5 mag dust 
correction will need to be tested in the future with follow-up near-IR spectroscopy. 

Steps (2), (3), and (4) are then performed by following the methods used in Ly et al. 
(2011), which we briefly review here. The effective volume is computed, accounting for the 
deviation of the filter profile from a perfect top-hat function. The maximum surveyed vol- 
ume is 3.83 X 10^ Mpc^ deg~^ for Ha emission-line fluxes above 2.7 x 10~^^ erg cm~^ s~^, 
and decreases to half of this volume at an Ha emission-line flux of 1.35 x 10"^'' erg cm~^ 
s~^. The total area covered by the stacked narrowband observations is 127 arcmin"^, and 
hence the resulting maximum volume is 1.36x10^ Mpc^. Next, Monte Carlo simulations are 
performed to estimate the completeness fraction at a given Ha luminosity. The simulations 
are used to compute the probability of detecting a source of a given continuum brightness 
and emission-line flux (i.e., NB210 excess), based upon the sensitivity of the narrow and 
broad-band images. We follow a maximum-likelihood approach where assumptions about 
the true distributions of the z ~ 2 star-forming population are varied to best reproduce the 
observed Ha LF and EW distribution. To model the intrinsic Ha EWs, log-normal Gaussian 
distributions are adopted, where the rest-frame median logarithmic EW and standard devia- 
tion are varied in 0.1 dex increments between 1.8 and 2.3; and 0.1 and 0.6, respectively. The 
Kg-hand luminosity function for 2 ~ 2 star-forming BzK galaxies, which is determined from 
data published in Ly et al. (2011b), is assumed. With the X^-band brightnesses and Ha 
EW's, 400,000 sources for each of 6x6 different EW model distribution are generated. Noise 
is then added to the sources, and the selection criteria described in Section 5.2 are applied 
to calculate the completeness as a function of Ha flux. The EW model that best reproduces 
both the observed Ha LF and EW distribution is described by a median log(EWrest) of 
2.1 and (T[log(EW)] =0.3 dex. The completeness of the sample is 99% at a flux of 2.2 x 
10"^'' ergs s~^ cm~^ and 50% at 1.2 x 10^^^ ergs s~^ cm~^. For comparison, the median 
log(EWrest) is 1.5 for local galaxies along the star-forming sequence (Lee et al. 2007). The 
characteristic specific SFR, which is traced by EW(Ha), is a factor of ~3 higher at z ~ 2 
(c.f., Noeske et al. 2007). 

Figure [TT] presents the Ha luminosity functions. The left panels show the LFs based 
on the dual narrowband excess selected sample (N=41), while the LFs in the right panels 
also include the NB210 excess galaxies with photo-z's between 1.8 and 2.6 (N=50). In every 
plot, there are two sets of points corresponding to (1) the raw Ha-|-[NII] luminosities, and 
(2) the fully-corrected luminosities, where the contribution of the [Nil] lines, 0.5 mag of dust 
attenuation, and survey incompleteness have all been taken into account. The data points 
for all four luminosity functions are given in Table [2l The bottom row is identical to the top 
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row, except that the results of Hayes et al. (2010; hereafter HIO) and Sobral et ah (2012; 
hereafter S12) are overplotted. 

Both the raw and fully-corrected LFs appear to rise steeply at the faint end. Single 
power-law fits to the faintest three, four or five bins for the four samples shown in Figure [TT] 
yield power law exponents in the range —1.60 < a < —1.93, with la uncertainties of ~ 0.3. 
These values are consistent with the HIO VLT/HAWK-I result of a = —1.77 ± 0.21, which 
was based on a photometric redshift selected 2.095 /im narrowband excess sample (N=55). 
Such results may imply that the star formation rate distribution function gradually steepens 
from the present day value of a = —1.2 zt 0.2 (e.g. Perez- Gonzalez et al. 2003), as has been 
suggested previously (e.g., Hopkins et al. 2000; Reddy & Steidel 2009; HIO; but see S12), 
though the uncertainties are still too large to conclude that evolution in a occurs. Robust 
measurements of this evolution will provide insight into hierarchical mass assembly through 
gas-rich mergers, and the regulation of star formation in low-mass halos. 

Of course, however, the value of a given by the power-law fit is an upper- limit. How 
closely this approximates the true value depends on how far the data points are from L*, 
the knee of the LP. Given the volume probed by our observations, and previous estimates 
of L* = io^3.07±o.22 gj,gg g-i (^fjxO; uncorrected for dust), it is clear that the current PourStar 
data do not provide useful constraints on the bright, exponential end. A complete determi- 
nation of the z = 2.2 Ha LP from the NewHa survey awaits the integration of wide-shallow 
data from NEWPIRM observations (C. Ly et al. 2012, in preparation). However, initial 
comparisons can be made with the recent study of S12, who combine VLT/HAWK-I and 
UKIRT/WPCAM imaging in the standard H2 narrowband filter at 2.121 /jm, and select their 
Ho; emitters with a mix of photometric redshift, BzK, and UBR selection criteria (N=45 
and 518, for the HAWK-I and WPCAM observations respectively). The bottom panels of 
Pigure [TT] show that, despite the patchwork of selection methods applied by S12, there is 
good agreement between the two studies at the faint end of the fully-corrected LP data- 
points. S12 report a = —1.57 ± 0.2 with their additional data covering higher luminosities, 
as they find a relatively low L* of io^^-^^='=''-^ ergs s~^ (for 0.5 mag attenuation and a fixed at 
-1.6). If the PourStar datapoints are combined with those from the S12 UKIRT/WPCAM 
imaging, the best fit Schechter parameters are L* = io^2.64±o.2i ^^^^ g-i^ log(0*)= -2.70±0.30 
and a = -1.58±0.40. 



7. Summary and Future Work 



In this paper, we have presented first results from near-infrared narrowband imaging 
observations with the new PourStar camera at the 6.5m Magellan telescope. We have shown 
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that our implementation of the dual narrowband technique for isolating deep intermediate 
redshift emission-line selected samples is highly effective, and identifies the z = 2.2 Ha 
emitters our NB210 excess sample with >80% completeness. The dataset is the deepest 
of its kind: the exquisite image quality (~0'.'55 FWHM) resulting from the combination 
of instrument, telescope and site results in 1.19 ^m and 2.10 /im imaging which probe 3a 
emission- line depths down to 1.0x10"^^ erg s~^ cm~^ and 1.2x10"^^ erg s~^ cm~^ respec- 
tively, and is competitive with recent HAWK-I narrowband observations at the 8.2m VLT 
(e.g., Hayes et al. 2010). Our z — 2.2 Ha sample extends the luminosities and SFRs probed 
by the well-studied UV-selected z ~ 2 galaxies of Erb et al. (2006c) by a factor of at least 
two. Optimization of the data reduction and analysis will further enable us to increase the 
depth of the sample by up to 40%. 

We compute constraints on the faint-end slope of the 2; = 2.2 Ha luminosity function 
using both the dual narrowband selected Ha emitters (N=41), and a combined sample of dual 
narrowband and photometric redshift selected galaxies (N=50). Fitting of a pure power-law 
gives an upper-limit on a of — 1.85±0.31, which is steeper than other recent estimates based 
on coarser selection techniques, but consistent within the large uncertainties that currently 
characterize such measurements. Combining our LF points with those at higher luminosities 
from other work, the slope decreases to a = —1.58 ±0.40. Additional data recently obtained 
from an FourStar pointing targeting the CANDELS area in the CDF-S will enlarge the 
sample, and tighten this constraint. 

These "narrow-deep" FourStar observations have been obtained as part of the larger 
NewHa Survey, which will combine the data with "wide-shallow" imaging through a similar 
narrowband filter pair with NEWFIRM at the KPNO/CTIO 4m telescopes, to carry out a 
statistical study of both luminous (but rare) and faint emission-line galaxies in the inter- 
mediate redshift universe. The dataset, in combination with the wealth of multiwavelength 
observations available in the extragalactic deep fields targeted, form the basis for follow-up 
studies on the physical properties (morphologies, dust, metallicities, recent star formation 
histories and stellar populations) of star-forming galaxies at a critical phase in their evolu- 
tion. The FourStar observations, in particular, enable study of lower mass members of the 
population, even probing SFRs down to LMC-type activities of a few tenths at 2; ~ 0.8. 

This project, including the purchase of the custom narrowband filters, was made possible 
through the support of Hubble and Carnegie Fellowships awarded to JCL. Special thanks are 
due to the FourStar instrument team and Carnegie Observatories, whose dedicated, long- 
term efforts to build an exceptional camera have enabled the science imagined for over a 
decade. 
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Table 1. FourStar Observations for a = 10''00'^32^4 6 = +02°16'"58" 



Filter 


Int. Time 


FWHM 


3-a Depth 




(hours) 


(") 


(AB, in 1'.'2 aperture) 


NB119 


9.65 


0.56 


25.1 


NB210 


8.39 


0.53 


24.2 


J2 


7.16 


0.61 


26.0 


J3 


8.36 


0.60 


25.7 


Ks 


5.82 


0.51 


25.0 



Table 2. Ha Luminosity Functions at z — 2.2 



log L 


$(L) 


N 


log L 


$(L) 


N 




Dual emitters (N==' 


41) 






Raw 






Corrected 




41.77 


6.838E-03 ± 3.419E-03 


4 








41.90 


5.861E-03 ± 1.566E-03 


14 


41.94 


8.096E-03 ± 2.649E-03 


9 


42.10 


4.058E-03 ± 1.223E-03 


11 


42.10 


5.027E-03 ± 1.451E-03 


12 


42.30 


2.582E-03 ± 9.760E-04 


7 


42.30 


4.429E-03 ± 1.278E-03 


12 


42.50 


1.476E-03 ± 7.378E-04 


4 


42.50 


2.214E-03 ± 9.038E-04 


6 


42.70 


3.689E-04 ± 3.689E-04 


1 


42.70 


7.378E-04 ± 5.217E-04 


2 


Dual emitters + photo-z selected (N=50) 




Raw 






Corrected 




41.74 


5.569E-03 ± 2.274E-03 


6 








41.90 


7.966E-03 ± 1.878E-03 


18 


41.93 


1.021E-02 ± 3.005E-03 


11 


42.10 


4.427E-03 ± 1.278E-03 


12 


42.10 


6.703E-03 ± 1.675E-03 


16 


42.30 


2.951E-03 ± 1.043E-03 


8 


42.30 


4.798E-03 ± 1.330E-03 


13 


42.50 


1.844E-03 ± 8.248E-04 


5 


42.50 


2.952E-03 ± 1.044E-03 


8 


42.70 


3.689E-04 ± 3.689E-04 


1 


42.70 


7.378E-04 ± 5.217E-04 


2 
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redshift of Ha redshift of Ha 

0.7 0.8 0.9 2.1 2.2 2.3 




wovelength (^im) 



Fig. 1. — Bandpasses of the NewHa narrowband filter pair (magenta) together with those 
for the Z-FOURGE medium-band filter set (blue), and the standard Ks broadband (black), 
which are used in FourStar. The atmospheric transmission is also plotted (gray). The 
expanded view shows the profiles for the narrowbands overplotted on the OH sky spectrum 
(gray). Similar sets of filters are also used in the NEWFIRM camera. The narrowband 
filters probe low OH airglow regions in the sky spectrum, and are sensitive to Ha at z = 0.8 
(near the beginning of the ten- fold decline in the cosmic SFR density) and at z = 2.2 (within 
the peak of the cosmic star formation history). Sky spectrum and atmospheric transmission 
data from Gemini Observatories. Filter transmission data (shown in the upper plots) from 
Barr Associates. FourStar/Magellan system throughputs (shown in the lower plot) from S.E. 
Persson et al. 2012, in preparation. 
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2.16 2.18 2.2 2.22 2.24 



redshift of galaxy 

Fig. 2. — The NcwHa narrowband filters are coupled such that the [Oil] emission of Ha 
emitters selected with NB210 observations can be captured by the NB119 filter. The redshift 
ranges over which the filters are sensitive to [Oil] and Ha in the same volume are illustrated 
in this Figure. Bandpass shifting over the field-of-view of the detector is not significant 
because of the positioning of the filters near the focal plane in the camera. The relative shift 
between the center and the corner of the field has been computed to be -6.IOA for NB119, 
and -5.49A for NB210, which are 6% and 2.6% of the bandpass widths of the two filters. The 
redshift range over which NB119 is sensitive to [Oil] completely encloses the range where 
NB210 will detect Ha. 
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Fig. 3. — Color-magnitude diagrams illustrating the selection of NB210 excess (left panel), 
and NB 2 10- detected, NB119 excess (right panel) sources, where the photometry has been 
measured in 1'.'2 diameter apertures. The upper horizontal line indicates values which are 
five times the scatter in the color for continuum objects with magnitudes between 19.0 and 
21.5. The curves indicate the values of the color excess that are significant at the 3-cr level. 
In both panels, the NB119 excess sources are marked with blue circles, and the NB210 
excess sources with red circles; hence the "dual emitters" are marked in both red and blue. 
Three stars (black crosses) are identified via inspection of SEDs constructed from broadband 
photometry. 30% of the NB210 excess sources are dual emitters. 
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Fig. 4. — Color-magnitude diagram illustrating the selection of a NB119 excess sample. 
Whereas the NB119 color-magnitude diagram in Figure [3] shows photometry only at locations 
where there are NB210 detections, in this Figure, all sources detected on the NB119 image 
are plotted. Again, the photometry has been measured in l'!2 diameter apertures; the upper 
cyan horizontal line indicates values which are five times the scatter in the color for continuum 
objects with magnitudes between 19.0 and 21.5; and the curve indicates values of the color 
excess that are significant at the 3-cr level. 
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Fig. 5. — Distribution of photometric redshifts for the NB119 excess selected sample. 
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Fig. 6. — Distribution of photometric redshifts for the NB210 excess selected sample. Dual 
NB210/NB119 excess emitters are shaded. 
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B-z' 

Fig. 7. — A BzK color-color diagram for all objects detected in the NB210 imaging. NB210 
excess selected objects and dual NB210/NB119 excess emitters are marked as indicated in 
the figure. Criteria used to distinguish different populations of galaxies as given by Daddi 
et al. (2004) are also shown. The dual excess emitters nearly all satisfy the criteria for 
classification as ^^sBzK" galaxies, indicating that they are star-forming galaxies at redshifts 
between 1.4 and 2.5, and demonstrating that the dual excess selection technique succeeds 
in identifying the z = 2.2 Ha emitters in the sample. Five dual excess emitters have photo- 
z's less than 0.4 (blue triangles), but their SEDs can also be reasonably fit with models at 
z = 2.2, as corroborated by their BzK colors. 
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Fig. 8. — Top: [OII]A3727 and Ha luminosities for local galaxies in the SDSS DR7 sample. 
The emission-line flux sensitivity limits of the NewHa FourStar narrowband imaging are 
illustrated by showing the SDSS galaxies whose Ha emission would be detectable in the 
NB210 data if they were at z = 2.2 (yellow and cyan regions), and those which would also 
have their [OH] detected in the NB119 data (cyan region). Median values computed in 
bins of L([011]) are plotted along with the la widths of the distributions. Bottom: The 
sample of NB210 excess emitters are overplotted (black crosses), and dual NB210/NB119 
excess emitters are indicated (blue open circles). Objects which have photo-z's which would 
make them Ha candidates are also marked (red open circles). The dual narrowband excess 
technique, as implemented here, is estimated to identify >80% of z = 2.2 Ha emitters in the 
sample. 
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Fig. 9. — Comparison of the z = 2.2 dual narrowband selected Ha emitters (blue circles) 
with the rest-UV selected z ~ 2 sample of Erb et al. (2006c) (black crosses). In both 
panels the plotted measurements for the Ha emitters include both Ha and [NII]AA6548,83, 
while those for the UV-selected sample only include Ha. A factor of two correction for 
spectroscopic sht losses has been accounted for in the hne measurements of the UV selected 
sample, as determined by Erb et al. (2006c). All quantities shown are as-observed, with no 
correction applied for dust attenuation. The left panel shows the rest-frame equivalent width, 
while the right panel shows the line luminosity, scaled by 7.9x10"^^ to give indicative star 
formation rates (Kennicutt 1998). Both are plotted as a function of the Kg AB magnitudes. 
On average, narrowband Ha sample extends the luminosities and SFRs probed for z ~2 
star-forming galaxies lower by a factor of at least two. 
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Fig. 10. — Tests of local empirical relations, which are frequently used to correct for dust 
reddening and the contribution of [Nil], using the measurements of Erb et al. (2006a,b,c). 
Left: The attenuation of Ha for the Erb UV selected sample is plotted as a function of 
observed SFR, and the curve shows the relation given by Hopkins et al. (2001). Erb et al. 
(2006b) provide E(B-V)'s derived by SED fitting, and these are scaled using the Calzetti 
et al. (2000) extinction law to give the attenuation in Ha, assuming that reddening in 
the gas is about twice that of the stars (black crosses) and that it the same as that of 
the stars (light blue triangles). In either case, the Hopkins prescription will over-predict 
the attenuation. Right: EW(Ha+[NII]A6583) is used to predict [NII]A6584/Ha for the 
Erb UV selected sample using the local relation of Villar et al. (2008). The results are 
plotted as a function of Ks magnitude, to enable comparison with the Erb et al. (2006a) 
stacked spectroscopic [Nil] /Ha measurements (red circles). Here the local relation yields 
reasonable predictions, albeit with large scatter. To correct for the contribution of [Nil] to 
the narrowband fluxes for the z = 2.2 Ha emitters, the relation shown by the dotted line is 
used to estimate the [NII]A6584/Ha ratio from the Kg magnitude. 
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Fig. 11. — Ha luminosity functions aX z = 2.2. The left panels show the LFs based on the 
dual narrowband excess selected sample (N=41), while the right panels show those based 
on the dual narrowband excess selected sample plus NB210 excess galaxies with photo- 
z's between 1.8 and 2.6 (N=50). In each plot, the raw observed LP (no corrections for 
[Nil], dust, or completeness) is indicated by the blue crosses, while the LF with all of these 
corrections applied is given by the black squares. Error bars for the former are shown slightly 
offset in L(Ha;) for clarity. The top and bottom rows are identical, except that results from 
other recent Ha narrowband studies are overlayed. LFs from Hayes et al. (2010) (magenta 
triangles), and from Sobral et al. (2012) (red curve: [Nil] correction applied]; magenta curve: 
[Nil] , dust and completeness corrections all applied) are shown. Dust attenuation corrections 
have been made consistent with those applied here. 



